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ABSTRACT 

Context. Be/X-ray binaries show outbursts with peak luminosities up to a few times 10^^ erg/s, during which they can he observed and 
studied in detail. Most (if not all) Be/X-ray binaries harbour accreting pulsars, whose X-ray spectra in many cases contain cyclotron 
resonant scattering features related to the magnetic field of the sources. Spectral variations as a function of luminosity and of the 
rotational phase of the neutron star are observed in many accreting pulsars. 

Aims. We explore X-ray spectral and timing properties of the Be/X-ray binary GX 304-1 during an outburst episode. Specifically, 
we investigate the behavior of the cyclotron resonant scattering feature, the continuum spectral parameters, the pulse period, and 
the energy- and luminosity-resolved pulse profiles. We combine the luminosity-resolved spectral and timing analysis to probe the 
accretion geometry and the beaming patterns of the rotating neutron star. 

Methods. We analyze the INTEGRAL data from the two JEM-X modules, ISGRI and SPI, covering the January-Eebruary 2012 
outburst, divided in six observations. We obtain pulse profiles in two energy bands, phase-averaged and phase-resolved spectra for 
each observation. 

Results. We confirm the positive luminosity-dependence of the cyclotron line energy in GX 304-1, and report a dependence of the 
photon index on luminosity. Using a pulse-phase connection technique, we find a pulse period solution valid for the entire outburst. 
Our pulse-phase resolved analysis shows, that the centroid energy of the cyclotron line is varying only slightly with pulse phase, while 
other spectral parameters show more pronounced variations. Our results are consistent with a scenario in which, as the pulsar rotates, 
we are exploring only a small portion of its beam pattern. 
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1. Introduction 


GX 304-1 is a Be/X-ray binary (BeXRB) system discovered as 
an X-ray source in 1967 during a balloon observation (Lewin 
St al. 1968a,b). Subsequently, the source was established to be 
an X-ray pulsar with a pulse period of ~ 272 s (McClintock 
St al. 1977). A study of the recurrent outburst activity revealed 
a ~ 132.5 d periodicity, likely due to the system’s orbital period 
(Priedhorsky & Terrell 1983). The optical counterpart of the bi¬ 
nary is a B2 Vne star, whose distance has been measured to be 
2.4 + 0.5 kpc (Parkes et al. 1980). Since 1980, the source entered 
an X-ray off-state (Pietsch et al. 1986), showing no detectable 
emission for 28 years. The quiescence was interrupted in June 
2008, when INTEGRAL detected hard X-ray emission from the 
source (Manousakis et al. 2008). Since then, GX 304-1 lighted 
up repeatedly, becoming a periodically outbursting X-ray source. 
The period of the outbursts after 2009 is roughly the same as be¬ 
fore 1980, i.e., ~132.5d. The peak luminosities are <10^^ erg/s 
in the 5 - lOOkeV energy band. The origin of the X-ray emis¬ 
sion is believed to be accretion of matter from the circumstellar 
equatorial disk around the optical companion onto a magnetized 
neutron star. The strong magnetic field of the accretor (~10^^ G) 
channels the captured matter towards its magnetic poles where 
X-ray emission originates in an accretion structure. 


A Cyclotron Resonant Scattering feature (CRSE), or cy¬ 
clotron line, has been detected in the spectrum of GX 304-1 with 
a centroid energy of ~ 52 keV (Yamamoto et al. 2011). CRSEs 


are important features in the spectra of accreting pulsars. In a 
strong magnetic field, electron energies corresponding to their 
motion perpendicular to the magnetic field lines are quantized 
in Landau levels, causing resonant scattering of impinging pho¬ 
tons. The first such line ever was detected in data from a balloon 
observation of Her X-1 (Triimper et al. 1978). Nowadays, they 
have turned out to be rather common in accreting X-ray pul¬ 
sars, with ~ 20 objects being confirmed cyclotron line sources, 
with several objects showing multiple lines (up to four har¬ 
monics in 4U 0115-1-63, Santangelo et al. 1999). Reviews are 
given by e.g. Coburn et al. (2002); Staubert (2003); Heindl et al. 
(2004); Terada et al. (2007); Wilms (2012); Caballero & Wilms 
(2012). The energy of the fundamental line Ecyc is directly pro¬ 
portional to the magnetic field strength at the emission site, 
Ecyc ~ 11.6 X Bi2(1 + Zg) keV, where B \2 is the magnetic field in 
units of 10'^ G, and Zg is the gravitational redshift. 

More recent observations have shown that the cyclotron line 
energy in GX 304-1 is positively correlated with the observed 
luminosity (Klochkov et al. 2012). Such a positive correlation 
was first observed in Her X-1 by Staubert et al. (2007) and is 
now observed also for GX 304-1 , Vela X-1 (Eiirst et al. 2014) 
and recently also for A 0535-1-26 (Sartore et al. 2015). 

The opposite correlation, a negative dependence between 
the cyclotron line energy and the luminosity, was actually de¬ 
tected earlier in high luminosity transient sources (4U 01 15 h- 63, 
Cep X-4, and V 0332-H53, Mihara et al. 1998). During the 
2004/2005 outburst of V 0332H-53, a clear anti-correlation of 
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Fig. 1. The SwiftfBKI light curve of GX 304-1 during the January - 
February 2012 outburst. The horizontal bars mark the INTEGRAL ob¬ 
servations. The numbers above the bars indicate the last digit ’i’ of the 
respective INTEGRAL revolution number - 113i. 


the line position with X-ray flux was observed (Tsygankov 
et al. 2006). A similar behavior was observed in outbursts of 
4U 0115+63 (Nakajima et al. 2006; Tsygankov et al. 2007). The 
presence of the correlation in this source, however, is shown to 
depend on the used spectral model (Muller et al. 2013; Iyer et al. 
2015). 

It has also been found that the power law spectral in¬ 
dex r (the absolute value) shows the opposite correlation with 
flux compared to the cyclotron line energy, both on shorter 
(Klochkov et al. 2011) and longer (Fiirst et al. 2014) time scales. 

A model-independent way to probe the correlation of the 
spectral hardness on flux is through a hardness-intensity diagram 
as shown by Reig & Nespoli (2013) for BeXRBs, which follow 
certain patterns in the hardness-intensity diagram and undergo 
state transitions if they reach a luminosity above a certain crit¬ 
ical value. For the physical interpretation of these correlations, 
that is the association with sub- and super-Eddington accretion 
regimes, see Becker et al. (2012) and Poutanen et al. (2013). 

Recently, GX 304-1 has passed through a new period of qui¬ 
escence, showing almost no X-ray activity between 2013 and 
2015. The data taken during the new active episode allowed Su- 
gizaki et al. (2015) to derive orbital elements for this binary sys¬ 
tem. Only recently has the source resumed its X-ray activity. 

Here we present the results of the timing and spectral stud¬ 
ies of GX 304-1 which was observed by INTEGRAL during an 
outburst in January-Eebruary 2012. Our work is primarily fo¬ 
cused on the timing analysis, which allowed to phase-connect 
the observations throughout the entire outburst, and on the pulse 
phase-resolved spectroscopy, with the goal to probe the emission 
geometry of the rotating NS at different viewing angles. The X- 
ray continuum of the source shows variability with pulse phase, 
with some features present only at particular pulse-phase inter¬ 
vals. Changes of the spectral parameters with pulse phase are 
quite common among accreting X-ray pulsars and are generally 
attributed to a change in the viewing angle of the emitting re¬ 
gion (see, e.g., Klochkov et al. 2008a and Vasco et al. 2013). 
Located relatively nearby and having a prominent cyclotron line 
feature, GX 304-1 is well suited for such studies. We confirm the 
results of Klochkov et al. (2012) using the latest version of the 
INTEGRAL analysis software and calibration. We investigate 


Table 1. INTEGRAL observations of the GX 304-1 outburst. 


Rev# 

Obs. ID 

Mid. MJD 

Exposure [ks] 

JEM-X ISGRI SPI 

1131 

09400230006 

55944.0 

64.6 

42.7 

68.6 

1132 

09400230007 

55947.0 

42.4 

31.9 

36.6 

1133 

09400230008 

55950.0 

- 

- 

10.7 

1134 

09400230009 

55952.8 

7.3 

25.4 

37.8 

1135 

09400230010 

55955.7 

- 

6.7 

25.1 

1136 

09400230011 

55958.7 

36.9 

28.1 

32.9 

1137 

09400230012 

55962.0 

78.1 

59.7 

78.4 

1138 

09400230013 

55965.0 

60.7 

45.2 

52.3 


the continuum variations (photon index and hardness) with lu¬ 
minosity. We obtain a pulse period solution that is valid through¬ 
out the observed outburst. The measured pulse periods are used 
to study pulse profiles in two energy bands and phase-resolved 
broadband (5 - lOOkeV) spectra. 


2. Observations and data 

The data used in this work are the same as those used by 
Klochkov et al. (2012). The observation log is given in Ta¬ 
ble 1. The analyzed outburst started on the 8th of January, 2012 
(Yamamoto et al. 2012). The International Gamma-ray Astro¬ 
physics Laboratory (INTEGRAL, Winkler et al. 2003) observed 
the event starting around MJD 55943.5, when the source flux in 
the 20 - 80 keV energy band was ~ 250 mCrab, until MJD ~ 
55965.5, when the flux dropped to ~ lOOmCrab. The peak of 
the outburst is reached aroun MJD 55953, with a flux exceeding 
1 Crab. INTEGRAL performed a total of eight observations, one 
per each satellite orbit (i.e., one every about three days), with 
a typical exposure of a few tens of kiloseconds per each obser¬ 
vation. The three high-energy instruments onboard INTEGRAL 
allow observations in broad enery ranges: the two Joint European 
X-Ray Monitor units (JEM-X, Lund et al. 2003) are sensitive in 
the range 3-35 keV; the Imager on Board the INTEGRAL Satel¬ 
lite (IBIS, Ubertini et al. 2003) is sensitive from ~ 20keV to a 
few MeV; and the Spectrometer onboard INTEGRAL (SPI, Ve- 
drenne et al. 2003) is sensitive in roughly the same energy range 
as IBIS. The duration of INTEGRAL observations are shown in 
Eig. 1, compared to the Swift/BAfT light curve of the outburst'. 
The actual observations are shortened due to solar activity, lead¬ 
ing to a reduction of the exposure times. 

Eor our analysis, we decided to use only data from those ob¬ 
servations where all of the three instruments were active, to en¬ 
sure a broad band coverage and a homogeneous analysis of all 
observations. Therefore, the revolutions 1133 and 1135 are not 
taken into account. Of the IBIS instrument, we only used the data 
from the Integral Soft Gamma-Ray Imager detector (ISGRI, Le¬ 
brun et al. 2003), which is sensitive in the 20 - 300 keV energy 
range. However, according to the recommendation of the instru¬ 
ment team, for the spectral analysis we restricted the duty energy 
band of each instrument to a “safer” range which is believed to 
be free of strong systematics: 5 - 30keV for the two JEM-X 
modules, 22 - 100 keV for ISGRI and 25 - 100 keV for SPI. 
The upper limit of lOOkeV for ISGRI and SPI data is due to the 
reduced photon flux of GX 304-1 above this energy. Data reduc¬ 
tion was performed with the version 10 of the Offline Science 


' http://swift.gsfc.nasa.gov/results/transients/weak/GX304-1/ 
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Table 2. Global timing solution obtained from phase-connection of IBIS/ISGRI (18 - 80keV) pulse-profiles and spectral parameters throughout 
the outburst obtained with OSA 10 plus additional gain correction (see text). Fluxes and luminosities are calculated in the range 3 - 80keV. 


Rev. ID 

Reference values 

1131 

1132 

1134 

1136 

1137 

1138 

/>[sl‘ 

274.9817 

275.1313 

275.0772 

274.9763 

274.8714 

274.8134 

274.7597 

T [MJD] 

55952.4592 

55943.98022 

55947.04829 

55952.76433 

55958.70659 

55961.99595 

55965.03900 

Fhixl{\0~^ erg s~^ cm~^) 

- 


1 ■57Q+0.004 

i.j/2'_qoo4 

9 41 9+0.006 
z..‘tij_0Q07 

1 991+0.004 
i.z,ji_ooo4 

0 609+”-^”- 

A 7QQ+0.002 

U.Z56_qpq2 

ergs-') 

- 

0.45 

0.95 

1.66 

0.85 

0.42 

0.19 

r 

- 

1.18/^ 

0.89/g;“ 

A Q9+0.09 
^•^■^-0.09 

-0.07 

1 77-10.O8 
-0.07 

1 56‘^^-”^ 

-0.03 

Etoiii [keV] 

- 

21-2/;:^ 

16.8/“ 


1 "7 "7-10.8 

^ ^-0.7 

21 3‘^^- 

Zl.O-i j 

30/? 

E,y,[keV] 

- 

55.0/g;/ 

56.6/“ 

59 . 3 /“ 

54 

C 4 4+0.8 
-^^•^-0.7 

50.6/11 

a-cyclkeV] 

- 

6.6/i;( 

0 c-bO.4 
°-^-0.4 

10 . 8 /“ 

”•■^-0.4 

^•-^-0.5 

c -2+0.8 
-^••^-0.7 

'/eye 

- 

11.9/li 

15 7+*-^ 

17“ 

10.5/1:? 

10.4/1;? 


■^red 

- 

0.9 

0.8 

1.3 

0.9 

1.0 

1.0 


* The reference values are the pulse period P valid at the time T in MJD. The uneertainties of the pulse periods given is of the order of 0.10 ms. 
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Pulse phase 

Fig. 2. The IBIS/ISGRI (18-80 keV) normalized pulse profiles obtained 
with the timing solution resulting from our phase-connection analysis. 
Each panel represents, from top to bottom, the successive INTEGRAL 
observations (from Rev# 1131 to 1138, shown in the top-left comer). In 
the top-right comer the luminosity of the source (3 - 80 keV) in units of 
10'^’ erg/s is reported. The thin straight line around phase 1.0 at the bot¬ 
tom of each panel (except Rev# 1138) approximately marks the phase 
range where the sharp features used to phase-connect the observations 
are found. 

Analysis (OSA) software^. Following the OSA user manual, we 
added a systematic error to the final count rates at a level of 1% 
for ISGRI, 3% for JEM-X, and 0.5% for SPI. In Section 4 also 
the previous version of the software (OSA 9) is used to reduce 

^ http://www.isdc.unige.ch/integral/analysis 
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Pulse phase 

Fig. 3. JEM-X 1 normalized pulse profiles in the 3 - 20keV energy 
range. Panels follow the same Rev# and luminosity stage as in in Fig. 2. 

the same data in order to compare the results and to explore the 
changes in the calibration, especially in the ISGRI energy scale, 
between the two versions. Details of the comparison between 
OSA 9 and OSA 10 are given in Appendix A. 

3. Timing anaiysis 

For the timing analysis, we extracted light curves with a time 
resolution of 3 s from JEM-X 1 (in the range 3 - 20keV) and 
IBIS/ISGRI (in the range 18 - 80keV). All time stamps were 
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55945 55950 55955 55960 55965 55970 

Time [MJD] 

Fig. 4. Pulse period as a function of time obtained by phase-connection. 
The dashed line represents a linear fit to the data points yielding a pulse 
period derivative of P = (-2.04 ± 0.01) 10~^. The period error bars are 
smaller than the symbols. 



corrected such that they refer to arrival times at the solar system 
barycenter. For each INTEGRAL observation the pulse period 
was determined by epoch-folding (Leahy 1987). Using these 
pulse periods, we generated ISGRI pulse profiles in three energy 
ranges; 20 - 40 keV, 40 - 60 keV, and 18-80 keV. The shapes of 
the profiles do not vary much within these three energy bands. 
Thus, to maximize the photon statistics, we chose the 18-80keV 
range for further analysis. At the time of preparation, correction 
for the orbital motion in the binary was not possible since the 
orbital parameters of the binary system were unknown. An or¬ 
bital solution for GX 304-1 has been furnished in the meanwhile 
by Sugizaki et al. (2015). In our analysis, however, phase con¬ 
nection has been succesfully achieved without taking the binary 
motion into account. In order to refine the pulse periods neces¬ 
sary to obtain the pulse profiles, we applied the phase-connection 
techniques (Deeter et al. 1981) to the 18-80 keV pulse profiles, 
as described in the following. 

Lor each INTEGRAL observation, a pulse profile was con¬ 
structed and used to determine an absolute reference time (in 
MJD), the arrival time, of the first pulse in the interval by mak¬ 
ing use of well defined ’sharp features’ in the pulse profiles, such 
as - see Lig. 2 - the peak near pulse phase zero and the minimum 
around pulse phase 0.1 or the time at the flux level centered be¬ 
tween the maximum and the minimum. 

If only the first derivative P of the pulse period P is taken 
into account, the expected arrival time f„ of the n-th pulse is 

t„^to+nPQ + ^n^PP. ( 1 ) 

Here, Pq is the pulse period at the reference time fo, n is the pulse 
sequence number, and P is the first period derivative. The pulse 
period at the time t can be determined as 

P(t) ^ Po + (t - to)P. (2) 

A consistent determination of pulse arrival times is formally 
only possible if the pulse shape does not change with time. 

In our observations, however, the pulse shape varies substan¬ 
tially between the intervals. Still, the main features of the pulse 
profile can clearly be recognized in all observations. The uncer¬ 
tainty of the corresponding pulse arrival times is given by the 


Lx/10^’’ [erg s i] 

0.1 1.0 



Fig. 5. Cyclotron line centroid energy E^yc as a function of the logarithm 
of flux in the 3-80 keV range. The dotted line is the result of a linear fit 
to the Ecyc - logio(Flux) data. The error bars indicate Icr-uncertainties 
(the flux uncertainties are smaller than the symbol size). The top x- 
axis shows the corresponding isotropic source luminosity assuming a 
distance of 2.4 kpc. 


Lx/10^^[ergs 

0.1 1.0 



Fig. 6. Variation of the photon index F versus the observed flux. The 
dashed line represents the linear fit to the data points. The luminosity 
scale on the top x-axis is as in Fig. 5. 


available statistics in the profile for the given shape and sharp¬ 
ness of the feature, and in our case corresponds to about 0.01 in 
phase. Using the determined pulse arrival times, the phase con¬ 
nection between the observations can be established. The pulse 
period solution (fo, Pq, P) is found from the fit of the measured 
pulse arrival times by formula (1). The uncertainties of Pq and 
P are calculated from the;^'^ contours at Icr c.l. This method al¬ 
lows the pulse period evolution to be measured with substantially 
higher precision compared to epoch folding. 

The obtained pulse period solution is presented in Table 2 
and visualized in Fig. 4. A decrease of the observed pulse period 
through the outburst is seen. 

The initial pulse period is Pq - 274.9817 + 0.0001 at MJD 
55952.4592 and the period derivative is L = (-2.04 + 0.01)xl0“^ 
(constant over the time of observation). This solution results to 
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Fig. 7. Phase-averaged spectra of four observations (labelled above each 
spectrum), from top to bottom, with decreasing luminosity. JEM-Xl 
(black) and ISGRI (red) data and folded model are shown. Spectra have 
been scaled for better visualization. Both the £cyc-flux and the E-flux 
correlations are clearly visible. 

be inconsistent with that found by Sugizaki et al. (2015). How¬ 
ever, such a difference is expected due their different approach, 
which also takes into account the orbital doppler effects. With the 
obtained final timing solution we folded the ISGRI light curves 
to produce a refined set of pulse profiles shown in Fig. 2. We 
note that there is still a slight phase shift between the pulse pro¬ 
files, up to A0 ~ 0.1, as can be determined, e.g., by the phase of 
the sharp minimum after the largest peak. This is probably due 
to variations in the pulse profile shape, which is most evident in 
the last observation (see Figs. 2 and 3). The same timing solution 
has been used to fold JEM-Xl light curves between 3 - 20keV 
and to generate the JEM-X pulse profiles for each of the six ob¬ 
servations as shown in Eig. 3. 

4. Spectral analysis 

Eor the following spectral analysis, we used the version 12.7.1 
of the XSPEC software (Arnaud 1996). Eollowing Klochkov 
et al. (2012), we fitted the spectrum of GX 304-1 during the ob¬ 
served outburst with a high energy exponential rolloff model (the 
cutof fpl component in XSPEC); 

F{E) oc exp[£/£(3) 

(where E is the photon energy, F is the photon index and 
Efoid is the e-folding energy of exponential rolloff in keV) and a 
multiplicative Gaussian absorption line (the gabs component in 
XSPEC) to account for the CRSE: 

where Ecyc is the cyclotron line centroid energy, tj and cr 
are the line depth (while the optical depth at line center is 
rj! and the width of the line, respectively. The 6.4 keV 

Ee Ka emission line is modeled with an additive Gaussian, even 
if it does not significantly improve the fit. The photo-electric ab¬ 
sorption at low energies is modeled with the wabs component 



Fig. 8.^^ -contour plots of the folding energy Efou and the photon index 
r for four INTEGRAL observations. The contours correspond to;^^^ -t 1.0 
(the projections of this contour to the parameter axes corresponds to 
the 68%-uncertainty for one parameter of interest), and;^^" -l- 2.3 (68%- 
uncertainty for two parameters of interest). The respective luminosity 
(3 - 80keV) in units of 10^’ erg s“' is indicated. 


in XSPEC. Using data from both JEM-X modules, ISGRI and 
SPI, we extracted pulse-phase averaged spectra for each obser¬ 
vation and explored the behavior of the best-fit parameters over 
the outburst. The results are reported in Table 2, and confirm the 
positive correlation of Ecyc with flux in GX 304-1. The obtained 
values of the cyclotron line energy show, however, a systematic 
offset with respect to those reported in Klochkov et al. (2012), 
which, in turn, is reflected in a higher slope of the fitted line, 
i.e., 6.51 + 1.21 keV/logio(erg/cmVs). The offset results from the 
new ISGRI calibration implemented in OSA 10. Eor technical 
details about this offset we refer the reader to Appendix A. The 
variation of Ecyc with flux is shown in Eig. 5. 

As mentioned in the Introduction, in addition to the Ecyc-flux 
correlation, some X-ray pulsars exhibit a correlation between the 
spectral continuum hardness and flux. Using the power law in¬ 
dex F (which describes the spectral hardness below the folding 
energy), an anti-correlation with flux is observed, i.e., the ab¬ 
solute value of the index decreases (the spectra become harder) 
with increasing flux. This is also true for GX 304-1 for which we 
find a negative correlation between the photon index and the flux 
as shown in Fig. 6. A selection of spectra is also shown in Fig. 7, 
where the spectral variation with luminosity is clear. 

For both correlations, Ecyc/ logio(flux) and F/ logio(flux), we 
have determined the Pearson’s correlation coefficients p and the 
corresponding one-sided probabilities p of obtaining the corre¬ 
lations by chance: pcyc = 0.92, pcyc ~ 0.01, and pr = -0.95, 
Py ~ 0.004, respectively. 

The folding energy Efoid also shows a negative correlation 
with the observed flux, similar to that of the observed photon 
index. To verify whether the anti-correlation between the fold¬ 
ing energy and the photon index is artificial (model-driven), we 
produced contour plots of the two parameters. The correlation 
between the folding energy and the photon index can be seen 
in Fig. 8. Altough the contours indicate some intrinsic coupling 
between Efo\d and F, the confidence regions corresponding to dif¬ 
ferent observations are well detached. Therefore, there must be a 
physical correlation between the two parameters. In this p\ot,x^- 
contours for four observations are also indicated at different lu- 
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minosities, showing an anti-correlation between the folding en¬ 
ergy and the observed luminosity. 

5. Phase-resolved spectral analysis 

As the neutron star rotates, at each time (or phase) we observe 
the emitting regions of the neutron star at a certain viewing an¬ 
gle. To explore the corresponding spectral variability with pulse 
phase, we have extracted the X-ray spectra of GX 304-1 in five 
equally spaced pulse-phase intervals (a finer binning would lead 
to larger uncertainties). For this, we have filtered the INTE¬ 
GRAL data with the Good Time Intervals (GTIs) corresponding 
to our phase bins. 

To obtain the phase-resolved spectra, we used the data from 
both JEM-X modules, ISGRI and SPI. To fit the phase-resolved 
spectra, we used the spectral model described in Section 4. The 
model provides an acceptable fit (the reduced chi-squared values 
are between 0.8 and 1.2) to the spectra from all but one 
phase bins (see Table 3). The spectra of the 0.0-0.2 phase bins of 
all revolutions, except the one from Rev# 1138, however, show 
large residuals with;^^^^^ > 1.5. Strong positive residuals appear 
around 35 keV. 

These spectra can be modeled with an additional Gaussian 
emission component (see Eig. 9). However, to enable a mean¬ 
ingful comparison of the spectral parameters in all phase bins, 
in Table 3 we report the results of the phase-resolved analysis 
obtained using the same model as in Section 4 (that is, with¬ 
out the additional Gaussian emission component to model the 
residuals around 35 keV). A comparison of the models with and 
without the inclusion of the additional Gaussian line is shown 
in Table 4. Alternatively, instead of the Gaussian emission line 
to model the bump, a Gaussian absorption line at lower energies 
can be used to fit the residuals around 35 keV, resulting in an 
equally good fit. In this case, the line centroid energy is about 
20keV. As it can be seen in Table 3, the shape of the continuum 
changes over the pulse phase. At the same time, the cyclotron 
line energy Ecyc shows only slight variation with pulse phase. 
However, the results in Table 3 show large uncertainties and the 
variation of many parameters with the pulse phase can not be 
firmly established. Indeed, even if the nominal value of a given 
parameter shows some variation, it is generally consistent with a 
constant value. Eor a deeper investigation of the spectrum vari¬ 
ability with pulse phase, we tried to improve the quality of the 
phase-resolved spectra by stacking together the observations at 
different fluxes. We therefore have extracted X-ray spectra in 12 
pulse-phase intervals. The number of phase bins has been cho¬ 
sen to ensure a good statistics in each interval, while the width 
of the bin has been kept equal to 0.075. In two cases, larger bin 
sizes are used to obtain a meaningful fit. Using this approach we 
loose information on the luminosity dependence of the spectra, 
and possibly introduce some systematic effects due to the stack¬ 
ing of observations with different spectral shapes. However, it 
results in an improved statistics, and allows us to study the flux- 
averaged dependence of the spectrum on pulse phase. 

In this analysis we excluded the data of Rev# 1138, which 
suffers of much lower photon statistics. The used model for spec¬ 
tral fitting is the same described in Section 4, which provides 
a good fit (Xjsd between 0.9 and 1.2) for all phase bins except 
three; 0-0.075, 0.075-0.225 and 0.675-0.75. Similarly to the 
0.0-0.2 phase-bins spectra of the individual observations, these 
phase bins require an additional component to get an acceptable 

which can be modeled by a Gaussian emission line around 
35 keV. Also here, the residuals in these three phase bins can al¬ 
ternatively be modeled by a Gaussian absorption line at lower 



10 20 50 

Energy (keV) 


Fig. 9. Top: The spectrum of the phase bin 0.0-0.2 from revolution 1132 
fitted with a power law/rolloff model including the cyclotron line and an 
additional Gaussian component to account for a bump around 35 keV 
(see text). Both JEM-X (red and green), ISGRI (black) and SPI (blue) 
data are used. Middle: The residuals for a fit with the model without 
the additional Gaussian component. Large residuals appear around the 
cut-off energy, fiofto/n.-The residuals after the inclusion of the Gaussian 
emission component to model the bump. 

energies, i.e., around 17 keV. Due to the presence of a very sim¬ 
ilar feature in the phase-resolved spectra of individual observa¬ 
tions, we consider it to be unlikely that such a feature is a result 
of systematic effects. However, to ensure a meaningful compar¬ 
ison of the results, we used the same model for all phase bins, 
i.e., the model without the additional Gaussian component. In 
order to compare the phase-resolved results with the pulse pro¬ 
file we also produced average pulse profiles in the ISGRI and 
JEM-X energy bands (18-80 and 3 - 20keV, respectively), ex¬ 
cluding the data from Rev# 1138. In Eig. 10 we show the best-fit 
results of the stacked phase-resolved spectroscopy, together with 
the ISGRI and JEM-X average pulse profiles. The Ecyc does not 
show significant variation with pulse phase (except for the phase 
bin 0.525-0.6), while the photon index E and the Efou confirm 
a variation which, however, is less than a factor of ~ 2. 

6. Discussion 

6.1. Timing results 

We have explored the behaviour of the pulse period in 
GX 304-1 during its January-Eebruary 2012 outburst using a 
pulse-phase connection technique. The measured period is used 
to construct pulse profiles in two energy bands: 3-20 keV (JEM- 
X), and 18 - 80keV (ISGRI), as shown in Eig. 2 and Eig. 3. An 
apparent decrease of the pulse period over the outburst is clearly 
seen. The measured period derivative is P = (-2.04 + 0.01) 10“^ 
(Eig. 4 and Table 2). This trend is most probably due to a 
combination of the Doppler effect caused by the orbital motion 
of the neutron star and an intrinsic spin-up caused by the transfer 
of angular momentum from the accretion disk to the neutron star. 

Considering the JEM-X (3 - 20keV) and the ISGRI (18 - 
80keV) energy bands, the pulse profiles of GX 304-1 are 
strongly energy dependent (except in the very end of the out¬ 
burst, i.e.. Rev# 1138), similar to the majority of accreting pul¬ 
sars. Thereafter, these two bands only will be considered to com- 
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Table 3. Best-fit parameters of the pulse-phase resolved spectra. Fluxes are calculated in the range 3-80 keV. 


Rev# 




1131 





1132 



Phase-bin 


0.0-0.2* 

0.2 - 0.4 

0.4-0.6 

0.6-0.8 

0.8- 1.0 

0.0-0.2* 

0.2-0.4 

0.4-0.6 

0.6-0.8 

0.8- 1.0 

r 



u.vu_qo5 

1 14+0.12 
^•^^-0.05 

1 74+0.16 
^"^^-0.17 


1 704-0.08 

i-^^-0.04 

0.875);“ 

1 4C+0.8I 

-0.21 

u.oo-o 15 

0.83!);1) 

£fold [keV] 


49!« 

19!| 
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21+) 

285) 

2751 

185? 

235) 

1351 

1551 

[keV] 
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565) 

565| 

56 4+0-9 
-0.8 

6O5) 

625) 

525) 

5751 

o-cyc [keV] 


4 3+^^ '^ 
^••^-0.6 

7 

'••‘-0.9 

° -1.2 
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^cyc 


10 
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26“^^ 
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' '”-0.015 
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0 621'^^ °^- 
u.uzi_0Qi3 

0.1375“)) 

1 419+0.021 
k-+k7-oo20 

1 794+0.022 
i.j7‘t-oo22 

^•^”‘“-0.019 

1 497+0.019 

^red 


1.8 

1.0 

0.8 

0.9 

0.8 

2.4 

1.0 

0.9 

0.8 

0.8 

Rev# 
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Phase-bin 


0.0-0.2* 

0.2-0.4 

0.4-0.6 

0.6-0.8 

0.8- 1.0 

0.0-0.2* 

0.2-0.4 

0.4-0.6 

0.6-0.8 

0.8- 1.0 

r 


0 59+0.06 

0 39‘^o.i2 
^•■^^-0.15 

0.885g:" 

”•^^-0.18 
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-0.18 

-0.03 

0.805);)) 

1 07+0-1* 
^•”'-0.18 
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£'foid [keV] 


1 7 74-O.6 
^ ^•'^-0.6 
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13+1 
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335? 
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185) 
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1851 

Ecyc [keV] 


58+1 

57+1 

625| 
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585| 
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o-cyc [keV] 
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4.85“ 
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z.uzu_QQ49 
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1 779+0019 
-0.021 

^•^”^-0.024 

1 1 co+0.020 

i.lJO_QQjg 
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1.5 

1.1 

1.0 

1.1 

1.2 

1.9 

0.9 

0.7 

0.9 

0.8 

Rev# 
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1138 



Phase-bin 


0.0-0.2* 

0.2 - 0.4 

0.4-0.6 

0.6-0.8 

0.8- 1.0 

0.0-0.2 
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0.8- 1.0 

r 


-0.04 
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-0.05 
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-0.14 

I 54+01'^ 

-0.15 

2 Q7+0.23 
-0,29 

1-82!);)) 
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-0.07 
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-0.12 
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185J 
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225) 
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-12 
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•7 1-4 
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Ecyc [keV] 
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o-cyc [keV] 


4 <5+0.6 
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J./-19 
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J -^-4 9 

-1.7 

3 4+°-^ 
■^•^-0.8 
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°-^-1.8 

Ucyc 


^•”- 1.4 
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- 4.3 
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-0.010 
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U.Z.UU_QQ,g 
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u.jj‘t_Qgj3 
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1.9 

1.0 

1.1 

0.9 

0.9 

0.9 

0.9 

1.0 
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Notes: The starred phase-bins mark those spectra that need an additional gaussian emission line around 35 keV to get an accettable lit. 


Table 4. Best-fit parameters of the first phase-bin spectra using a model with and without additional gaussian emission line. 



Rev# 1131 

Rev# 1132 

Rev# 1134 

Rev# 1136 

Rev# 1137 

Parameter 

No Blimp 

Bump 

No Bump 

Bump 

No Bump 

Bump 

No Bump 

Bump 

No Bump 

Bump 

r 

1 70+00-1 
^•'”-0.04 

2 15+0.21 

“^- 0.20 

1 79 + 0.08 
^•'^^-0.04 

2 07+'^-^’^ 
^•”'- 0.54 
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”••'^-0.06 
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1 56^®°^ 
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A Oir+0.16 
2.26_„,5 
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17 
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56 
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-'^•^- 1.0 

^”•'- 0.6 
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-'^••‘- 0.7 

49 . 95 ?) 

(Tcyc [keV] 

A 7+0.7 
^•■^-0.6 

7 9 + 0.8 
•^•^- 1.0 

c 7+0.6 
-0.5 

^•”-0.5 

6.8!);) 

165) 

4 - 8 !);) 

A 1 q+0.69 
”-^°-0.09 

A ^+0.6 
^•-'-0.7 

5-8!?:? 

7/cyc 

10.5!!;) 

AQ+181 

68-36 

21‘*'^ 

205 )° 

115 ) 

245)= 

11-0!);) 

unconstrained 

9 A+1-5 
^•”-1.4 

90+263 

^”-48 

^bump [keV] 

- 

365) 

- 

255) 

- 

35!) 

- 

74 7 + 0.7 
•^^"^- 0.7 

- 

38!) 

^bump [keV] 

- 

125 ) 

- 

145 ) 

- 

7.5!“ 

- 

”•^- 0.5 

- 

1051 

normburnp" 

- 

n nnA+0,002 

0.006_qqq2 

- 

”•”^^-0.004 

- 

U.UJ_QQ4 

- 

U.UWl J_Q QQQ2 

- 

U.UUJ_Q0QJ 

Xljd.o.f. 

1.8/136 

1.0/133 

2.4/136 

0.7/133 

1.5/136 

1.0/133 

1.9/136 

0.8/133 

1.9/136 

0.9/133 


Notes. XSPEC normalization of the bump component, in units of photon keV ‘ cm ^ s * at 1 keV. 


pare our pulse profiles results. The ISGRI pulse profiles pre¬ 
sented in Fig. 2 show two sharp peaks and a hint of a third peak 
in between (around phase 0.5) which is not present in all ob¬ 
servations. In the JEM-X pulse profiles, however, all three peaks 
are clearly visible in all observations, with the exception of Rev# 
1138. Comparing the pulse profiles in the two energy bands, it 
appears that the most prominent of the peak in the JEM-X pulse 
profiles (peaking around phase 0.5 in the top panel of Eig. 3) 
tends to disappear at higher energies. While there is no com¬ 
monly accepted interpretation of the pulse profile dependence 
on energy and luminosity, we note that a similar behavior is ob¬ 
served in Her X-1 (Klochkov et al. 2008b), 4U 01154-63 and 


V 03324-63 (Tsygankov et al. 2007). A qualitative interpretation 
of the pulse profile energy-dependence (based on the results of 
Mowlavi et al. 2006) has been proposed by Tsygankov et al. 
(2007) to explain the variation of V 03324-63 pulse profiles in 
a purely geometric fashion: if the axis of the magnetic dipole 
is at an angle to the rotation axis (a necessary condition to ob¬ 
serve pulsation), the neutron star itself can eclipse some of the 
lower parts of the accretion columns, which are thought to pri¬ 
marily emit the harder photons. Thus, compared to the pulse pro¬ 
file at softer energies, the harder one lacks of (at least) one peak. 
On the contrary, the softer photons are primarily from regions 
higher up in the accretion column, and therefore visible during 
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Pulse phase 


Fig. 10. The best fit parameters of the phase-resolved spectra using the 
five stacked INTEGRAL observations. The vertical dashed lines indi¬ 
cate the phase-bins. The solid blue curve is the JEM-X average pulse 
profile, while the green curve is the ISGRI average pulse profile. The 
vertical error bars are Icr uncertainties, while the horizontal bars repre¬ 
sent the phase bin width. 


most of the spin period. On the other hand, gravitational light 
deflection also plays a crucial role in shaping the observed pulse 
profile. Indeed, depending on to the geometry of the emitting re¬ 
gion and the observer’s line of sight, this effect could make the 
accretion columns visible from all directions (Kraus 2001; Kraus 
et al. 2003), thus affecting the pulse profile in a significant way. 
In addition, the pulse profiles luminosity dependence probably 
reflects a change in the physical condition of the accretion col¬ 
umn. Indeed, for sub-critical sources, the height of the accretion 
column is anti-correlated with the accretion rate, and this could 
be at least partially responsible for the luminosity dependence 
of the pulse profiles. However, this is only a qualitative scenario 
to account for the observed pulse profile changes. More accurate 
descriptions of the pulse profiles dependences need more com¬ 
plicated models in order to account for the magnetic field config¬ 
uration, the scattering cross section along the accretion column, 
the shape of the emission beam and its dependence on luminos¬ 
ity and energy, and other parameters. On the other hand, other 
accreting X-ray pulsars show the opposite behavior, where the 
low-energy pulse-profiles are absorbed, and the extra peak only 
appears at higher energies (see, i.e., Kreykenbohm et al. 2004, 
for the case of GX 301-2). 

Although the overall morphology of the profile remains the 
same in all observations (which permits phase connection of the 
observations, see Sec. 6.3), the relative intensities of the peaks 
slightly vary with luminosity. Furthermore, in the last stage of 
the outburst (Rev# 1138), the pulse profile resembles a single 
broad peak, and a weak energy dependence. A similar flux and 
energy dependence of the pulse profile is also observed by Deva- 
sia et al. (2011) during a GX 304-1 outburst in August 2010, 
when the source reached a peak flux of ~ TOOmCrab (that is, 
~30% less than the outburst analyzed in this work). However, 
the different pulse profiles of Rev# 1138 likely indicate a pos¬ 
sible change in the beaming pattern from an almost purely pen¬ 
cil beam at lowest luminosities, to a mixed pencil-fan beam as 


the luminosity increases (see Becker et al. 2012, and references 
therein). 

The pulsed fraction is only ~20% in both energy bands, and 
shows only slight variation along the outburst, uncorrelated with 
the flux. This is contrary to many other X-ray pulsars, where a 
pulsed fraction > 50% is observed, correlated with energy and 
flux (Lutovinov & Tsygankov 2008). This again may be under¬ 
stood in a qualitative manner following a geometrical interpreta¬ 
tion: it is possible that, as the star rotates, we are exploring only 
a small area of the beam pattern (i.e., the observed flux from one 
emission region as a function of the direction to the observer), 
which in turn leads to an observed low pulsed fraction. This sce¬ 
nario is supported by the results discussed in Section 6.3. 


6.2. Phase-averaged spectral results 

With the new INTEGRAL calibration (Caballero et al. 
2013) and analysis software (OSA 10), we confirm the 
luminosity-dependence of the cyclotron line energy Ecyc in 
GX 304-1 reported in Klochkov et al. (2012), with a better re¬ 
striction of the spectral parameters. The observed £cyc shows a 
systematic shift of ~5% with respect to the energies reported 
in the previous work. As a consequence, the correlation found in 
this work has a slightly higher slope than that found by Klochkov 
et al. (2012). As mentioned in Section 1, such dependency has 
already been found in several cyclotron line sources. The Ecyc- 
flux correlation is positive in Her X-1 (Staubert et al. 2007), 
Vela X-1 (Furst et al. 2014), A0535-H26 (Sartore et al. 2015), 
and some indications suggest also in GX 301-2 (La Barbera 
et al. 2005). Pulsars with a confirmed negative Lcyc-flux corre¬ 
lation are V 0332H-63 (Tsygankov et al. 2006), and 4U 0115-f63 
where the negative correlation was found using the NPEX model 
(Nakajima et al. 2006), but not when using a cutoff-power law 
model (Muller et al. 2013), or when the fundamental is modelled 
by two independent sets of CRSF lines at ~11 keV and ~15 keV 
(Iyer et al. 2015). Thus, GX 304-1 further enlarges the group of 
sources showing the positive correlation. 

Klochkov et al. (2012) interpret this correlation according to 
the model described by Staubert et al. (2007) and Becker et al. 
(2012). In these models, a critical luminosity Lc discriminates 
between two accretion regimes: one in which the height of the 
emitting column is governed by the ram pressure of the infalling 
material (for L < Lc, i.e., “sub-critical sources”), and the other 
where a radiation pressure dominated emitting column is formed 
(for L > Lc, i.e., “super-critical sources”). As the luminosity 
L varies, the X-ray emitting height also moves within the ac¬ 
cretion column, thus the local magnetic field changes, leading 
to a variable cyclotron line energy with luminosity. The criti¬ 
cal luminosity below which a positive L’cyc-luminosity corre¬ 
lation is expected, has been calculated by Becker et al. (2012) 
as a function of the magnetic field and other physical param¬ 
eters of the NS. For typical neutron star parameters, it results 
in Lc ~ 1.5 X 10^’B[2^‘^erg s^'. For GX 304-1, this leads to 
~ 8 X 10^^ erg s *. During the analyzed outburst, the source lu¬ 
minosity clearly remains below this value (see Figs. 5 and 6). 
Thus, the observed positive £cyc-flux correlation is consistent 
with the expected behavior for a sub-critical source according 
to the physical picture outlined in Becker et al. (2012). More 
recently, a new calculation of the Lc function has been proposed 
by Mushtukov et al. (2015), where a more detailed accretion sce¬ 
nario is taken into account. The observed positive Fcyc-flux re¬ 
lation in GX 304-1 is consistent with their model as well. 
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Fig. 11. Hardness/intensity diagram for GX 304-1 during the analyzed 
outburst. The soft color definition is given in the text. The errors are Icr 
uncertainties. The errors on the ISGRI rates are smaller than the symbol 
size. The dashed red curve corresponds to the rising part of the outburst, 
the dotted blue one - to the decay. 

We have also investigated the correlation between the spec¬ 
tral photon index F and luminosity. In many accreting X-ray pul¬ 
sars, r is luminosity dependent. It had been found that in accret¬ 
ing pulsars showing a positive iscyc-flux correlation, the correla¬ 
tion between F and the flux is negative and vice versa (Klochkov 
et al. 2011). The exact mechanism of the continuum/flux depen¬ 
dence is not yet clear, especially for super-critical pulsars (with 
L > Lc)- In sub-critical sources, including GX 304-1 , the ob¬ 
served hardening of the spectrum with flux (Fig. 6 ) can be quali¬ 
tatively understood assuming that the pressure and electron den¬ 
sity inside the emitting structure increase with increasing mass 
accretion rate and luminosity. Indeed, such an increase should 
lead to a higher optical depth t, thus to a larger Comptonization 
parameter y - TkTel(meC^). The spectrum formed by Comp¬ 
tonization of the seed photons from the base of the column would 
then become harder. 

Finally, we have investigated the dependence of the folding 
energy Efoid on the photon index F and, at the same time, on lu¬ 
minosity. These correlations can be seen in the contour plot in 
Fig. 8 . The plot shows that isfoid is positively correlated with F, 
and negatively correlated with flux (£foid increases at lower lumi¬ 
nosities). A similar behavior is reported for other two accreting 
pulsars: A 0535-H26 (Muller et al. 2013) and RX J0440.9-I-4431 
(Ferrigno et al. 2013). These correlations can be qualitatively un¬ 
derstood as follows. It is generally accepted that the powerlaw- 
cutoff spectrum of accreting pulsars is a result of Comptoniza¬ 
tion, i.e., scattering of soft photons off hot electrons supplied 
by the infalling plasma. As the luminosity increases towards L^, 
radiation field begins to affect the accretion flow and Compton 
cooling becomes more efficient. This might lead to the observed 
shift of Zifoid reflecting the electron temperature Tg towards lower 
energies while, at the same time, the Comptonization process 
makes the spectrum harder. 

The dependence of the spectral hardness on flux has recently 
been investigated for a number of accreting pulsars in BeXRBs 
by Reig &. Nespoli (2013). They found that the X-ray emission 
of the sources follows certain patterns in the hardness/intensity 
diagram. At lower fluxes, the emission first becomes harder with 
increasing flux. At a certain flux a “turn-over” occurs, so that at 


higher fluxes the emission again becomes softer (Fig. 6 in Reig 
& Nespoli 2013). The luminosities corresponding to the fluxes 
of the turn-over are in all cases close to the critical luminosity Lc 
discussed above. Thus, the turn-overs in the hardness/flux depen¬ 
dences may reflect a transition between the aforementioned sub- 
and super-critical regimes characterized by the opposite E^yc- 
flux correlations. To compare our results with those presented in 
Reig & Nespoli (2013), we constructed a hardness/intensity dia¬ 
gram for GX 304-1 (Fig. 11). We use the same definitions of the 
“soft-color” (SC) - the ratio between the fluxes in the 7 - 10 keV 
and 4 - 7keV ranges. As expected in a sub-critical source, 
the luminosity corresponding to the turn-over is not reached in 
GX 304-1 and the diagram shows only the part where the hard¬ 
ness increases with luminosity, i.e., the “horizontal branch” ac¬ 
cording to the terminology used in Reig & Nespoli (2013). 

6 . 3 . Phase-resolved spectral analysis 

We modelled the phase-resolved spectra with a rolloff power-law 
model and a multiplicative Gaussian absorption line to account 
for the CRSF The maximum value of the line energy, E^yc ~ 
62 keV, is reached at the peak of the outburst. Rev# 1134, in the 
pulse phase bin 0.4-0 .6 (see Table 3). Using a value of Zg = 0.3 
for the gravitational redshift (assuming canonical neutron star 
parameters), this corresponds to a magnetic field strength of 6.9x 
10'^ Gauss. This is among the strongest observed magnetic fields 
in an accreting X-ray pulsar. 

The centroid energy of the cyclotron line shows only slight 
variations with pulse phase, not exceeding ~16% (with respect 
to 0.2 wide phase bins). This is contrary to some other X-ray pul¬ 
sars (see, e.g., the cases of Cen X-3, Burderi et al. 2000; GX 301- 
2, Kreykenbohm et al. 2004; and Her X-1, Klochkov et al. 2008a; 
Vasco et al. 2011), where the energy of the cyclotron line shows 
up to 30% variation with pulse phase. GX 304-1 , however, is 
not the only accreting X-ray pulsar where the cyclotron line en¬ 
ergy varies only by < 20 % with pulse phase: see, e.g., the case 
of 4U 1538-52 (Coburn et al. 2002). 

We find that the standard rolloff/power law model does not 
lead to good spectral fits for the 0.0-0.2 phase bins. These spec¬ 
tra need an additional component to obtain an acceptable fit. We 
modelled this component with a Gaussian emission line around 
35 keV. Further analysis with more physically motivated models 
are needed to investigate the feature in more detail, and to un¬ 
derstand its origin. The extra component is reminiscent of the 
"lOkeV feature" observed in the spectra of many accreting X- 
ray pulsars (Coburn et al. 2002; Klochkov et al. 2008a; Muller 
et al. 2013), for which no accepted explanation of its physi¬ 
cal origin exists. However, the centroid energy of the feature 
in GX 304-1 is at a much higher energy, ~ 35keV. Interpreted 
within this context, the feature is likely due to the inadequacy 
of the underlying phenomenological model. Alternatively, we 
found that the same feature can be fitted by a Gaussian absorp¬ 
tion line at lower energies, instead of an emission line. In this 
case, the absorption line centroid energy is between 12 keV (with 
o- = 14keV, 77 = 37), and 25keV (with cr = 9keV, 77 = 8 ). To 
explore the possibility of the feature being the fundamental cy¬ 
clotron line, we calculated the harmonic ratio between the con¬ 
firmed cylotron line at ~55 keV, assumed to be the first harmonic, 
and the centroid energy of the absorption feature, and found that 
it significantly deviates from the canonical ratio of 2. For the 
feature at 12keV, the ratio results in an anomalously high value 
of ~4.6, while for the feature at 25 keV the ratio is halved to a 
value of ~2.3. Relatively high deviations from an integer value 
of the harmonic ratios are observed also in other sources, e.g.. 
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Vela X-1 (Fiirst et al. 2014, and references therein), where the 
ratio between the confirmed harmonic and fundamental line en¬ 
ergies spans in the range 2.1-2.7. Those authors invoke photons 
spawning from the harmonic lines as responsible to influence the 
shape of the fundamental line, shifting its measured centroid en¬ 
ergy from its real value. However, values of the harmonic ratio 
as high as 4.6 have never been observed for the first harmonic 
line, and seem to exclude the possibility of the absorption fea¬ 
ture being the fundamental cyclotron line. 

Thanks to the successful phase-connection (see Section 3), 
we have been able to identify phase bins at different stages of the 
outburst and, therefore, to perform phase-resolved spectroscopy 
in order to explore the spectral variability with pulse phase. In 
particular, we have been able to perform perform phase-resolved 
spectroscopy stacking together all observations (with the excep¬ 
tion of Rev# 1138), in order to improve the statistics and investi¬ 
gate the pulse phase variability of spectral parameters with better 
precision (12 phase bins). The results are shown in Fig. 10 for 
three parameters of interest: the cyclotron line centroid energy, 
the photon index and the folding energy. To compare the phase- 
resolved parameters of the stacked analysis, with the pulse pro¬ 
files, we produced averaged JEM-X and ISGRI pulse profiles, 
which are also shown in Fig. 10. The value of Fcyc is constant 
at ~ 56 keV (with one phase bin showing a ~10% higher value, 
however, not statistically significant). The photon index F and 
the folding energy Efoia show more prominent variations with 
pulse phase, however, less than a factor of ~2. A similar varia¬ 
tion of the continuum parameters is observed also in other ac¬ 
creting pulsars (see, e.g.. Her X-1, Vasco et al. 2013; GX 301-2, 
Suchy et al. 2012; Vela X-1, Maitra 8c Paul 2013), but a consis¬ 
tent theoretical scenario which supports such variations is still 
missing. 

Similarly to the phase-resolved analysis of the separated ob¬ 
servations, a few phase-bins of the stacked phase-resolved spec¬ 
tra need an additional component to give an acceptable We 
modelled this component by a Gaussian emission line around 
35 keV. An equally good fit is obtained with an absorption Gaus¬ 
sian line at lower energies, i.e., between 16 keV (with cr - 9 keV, 
7] = 17), and 21 keV (with cr = 13 keV, rj - 24). 

Since our timing analysis shows a low pulsed fraction which 
is flux- and energy-independent, and the stacked phase-resolved 
analysis shows a remarkably stable Ecyc, we are tempted to ex¬ 
plain these results within the context of the same scenario. In¬ 
deed, the steadyness of the cyclotron line energy along the pulse 
phase can be understood as follows. According to model calcu¬ 
lations, the cyclotron line parameters should strongly depend on 
the viewing angle and, therefore, on the pulse phase (see Schon- 
herr et al. 2007, and references therein). However, if the geo¬ 
metrical configuration is such that we are exploring only a small 
portion of the beam pattern, then the parameters result mostly 
insensitive to the pulsar rotation. Therefore, the viewing angle 
with which we observe the accretion column does not vary ap¬ 
preciably with the NS rotation and, as a consequence, the Ecyc 
results almost independent from the rotational phase. 


7. Summary and conclusions 

In this work, we present the results from spectral and tim¬ 
ing analysis of INTEGRAL data of the accreting pulsar 
GX 304-1 observed during an outburst in January-Eebruary 
2012 . 

We studied the behavior of the source spectrum as a func¬ 
tion of the luminosity analyzing the pulse phase-averaged and 


phase-resolved spectral parameters. Our results confirm the posi¬ 
tive Ecyc-flux correlation in GX 304-1 with the new INTEGRAL 
calibration, and show a negative correlation between the pho¬ 
ton index and the luminosity. The positive Ecyc-flux correlation 
is also observed in other accreting pulsars whose luminosity is 
below Lc, like Her X-1 and Vela X-1, while the negative photon 
index-luminosity correlation has also been observed in the pulse- 
to-pulse analysis ofHerX-1 and A 0535H-26. Then, GX 304-1 is 
one of the few sub-critical sources with a confirmed photon 
index-luminosity negative correlation. To test whether the be¬ 
havior of GX 304-1 below in a model-independent way, we 
constructed a hardness/intensity diagram and verified that the 
turn-over point, which reflects the critical luminosity, is never 
reached, and the diagram only shows the so-called horizon¬ 
tal branch. Eurthermore, a correlation between the folding en¬ 
ergy and the photon index has been found, with a relative anti¬ 
correlation between the folding energy and the observed lumi¬ 
nosity. These correlations seem to fit well in the context of an 
emergent spectrum shaped by both bulk and thermal Comp- 
tonization of blackbody seed photons radiating from the accre¬ 
tion mound (where thermal Comptonization becomes important 
when the luminosity increases towards L^). In this context, the 
observed spectral changes are a consequence of the increasing 
Compton cooling efficiency at higher luminosity. 

Thanks to the pulse profiles main features, and to their clear 
identification in all observations (except the last one), we have 
been able to phase-connect the pulse profiles, which allowed us 
to obtain a pulse period solution valid over the entire outburst 
with remarkable precision, and to identify phase bins at different 
stages of the outburst. We measured a period derivative seen over 
the full outburst, which is probably a combination of the Doppler 
shift due to the orbital motion of the neutron star and an intrinsic 
spin-up due to the angular momentum transfer by the accreted 
material. 

The observed pulse profiles compared in the JEM-X (3 - 
20keV) and ISGRI (18 - 80keV) energy bands show a strong 
energy and luminosity dependence during the outburst, except at 
the end of the outburst (Rev# 1138), where the pulse profiles are 
similar in the two energy bands and show a much less complex 
morphology. While the energy dependent pulse profiles can be 
justified by a geometrical effect due to the rotation of the neu¬ 
tron star, the luminosity dependence probably reflects changes 
in the physical details of the accretion process in the sub-critical 
regime. Moreover, we observe only a low pulsed fraction, inde¬ 
pendent of energy and luminosity. 

The phase resolved analysis shows that the cyclotron line is 
detectable at all rotational phase bins of the neutron star. The 
maximum reached centroid energy (~62 keV) enabled us to infer 
one of the strongest magnetic field strength on an accreting X- 
ray pulsar (while the highest belongs to GRO J1008-57 with a 
cyclotron line energy of ~80keV, Bellm et al. 2014). 

Eurthermore, we observe no significant variations of the cen¬ 
troid energy with the pulse phase. On the other hand, the depth 
and the width of the line, show more pronounced variations with 
pulse phase, up to a factor of ~3, with no correlation with the 
pulse profile. In the stacked phase-resolved analysis we have 
been able to show how the cyclotron line energy is essentially 
constant along the pulse phase. 

All these considerations lead us to a possible qualitative sce¬ 
nario where the geometrical configuration of the observed NS is 
such that we are able to explore only a small part of the beam 
pattern. 

We found that, contrary to the phase averaged spectra, a few 
phase-resolved spectra lead to large residuals when fitted by a 


Article number, page 10 of 12 


C. Malacaria et al.: INTEGRAL observations of GX 304-1 


standard rollofF/power-law model: they need an additive com¬ 
ponent that we modeled by a Gaussian emission line around 
~35 keV, but that can be equally modeled by an absorption line 
at ~20 keV. The additional Gaussian emission line heavily modi¬ 
fies the resulting best-fit parameters, often leading to unphysical 
values, most likely indicating that the feature is due to the inad¬ 
equacy of the underlying phenomenological model. 
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Appendix A: The OSA 10 data reduction 
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Fig. A.l. Cyclotron line centroid energy £cyc as a function of the log¬ 
arithm of flux in the 13 - 80keV range. The black dots are the results 
of Klochkov et al. (2012), obtained with OSA 9 and corrected for the 
ISGRI energy gain drift. The red triangles are the results of this work, 
obtained with OSA 10 and further gain corrected. The dashed and dotted 
lines are the results of linear fits to the Ecyc-logio(Flux) data for the two 
sets resulting from OSA 9 and OSA 10 data reduction, respectively. The 
error bars indicate Icr-uncertainties (the statistical flux uncertainties are 
smaller than the symbol size). The top x-axis shows the corresponding 
isotropic source luminosity assuming a distance of 2.4 kpc. 

Our spectral analysis has been performed with the OSA 10 
software which, at the time of writing, is the latest version of 
the INTEGRAL analysis software. OSA 10 is supplied with a 
new energy calibration of the ISGRI detector, with respect to the 
previous version (OSA 9). As a result, a more reliable photon en¬ 
ergy reconstruction is ensured. The details about the new energy 
reconstruction are reported in Caballero et al. (2013). 

To investigate the difference in energy calibration between 
OSA 9 and OSA 10, we extracted the data set used in this 
work with both the software versions. Following Klochkov et al. 
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(2012), after the OS A 9 data reduction, we performed an ad¬ 
ditional gain correction based on the background spectral line 
complex of Tungsten (W, nominally at 58.8297 keV). After the 
OS A 10 data reduction, which includes a more accurate gain cor¬ 
rection, a small deviation of the inferred energy from the nom¬ 
inal Tungsten line was still noticeable. The deviation amounted 
up to h- 0.5 keV in the Science Windows (ScWs, i.e., INTEGRAL 
pointings) at the beginning of each INTEGRAL revolution. Ac¬ 
cording to the ISGRI instrument team (priv. comm.), the devi¬ 
ation can be attributed to possible orbital drifts in the detector 
response, which are responsible for systematic scatter when rev¬ 
olutions start. However, even if orbital variations are not cur¬ 
rently corrected, they are mitigated when averaged over a num¬ 
ber of revolutions. Nonetheless, to eliminate this deviation, we 
performed an additional gain correction, in order to ensure the 
stability of the Tungsten reference line. After the correction, the 
inferred Tungsten line energy in all ScWs is found to be consis¬ 
tent with its nominal value within 3cr. 


To illustrate the effects of the different energy calibrations, 
we plotted the Ecj^c-flux relation for the spectra obtained with 
the two software versions (Eig. A.l). We note that the energy 
range indicated in Klochkov et al. (2012) for the calculation of 
the flux is erroneously reported as 4 - 80 keV, while the correct 
range is 13 - 80keV (where the lower limit is constrained by 
the response matrix energy range of the ISGRI instrument). The 
two different sets of results clearly show a systematic shift of 
the cyclotron line centroid energy, with a deviation (proportional 
to the measured energy) up to ~2 keV. This is the result of the 
new energy calibration implemented in OS A 10 plus a minor 
contribution due to the additional gain correction performed after 
the OS A 10 data reduction. In our scientific analysis, we used the 
data reduced with OS A 10 with the additional gain correction. 
This ensures the most accurate measurements of the cyclotron 
line centroid energy. 
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